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1 Introduction 

In the Standard Model (SM) of particle physics, the Yukawa couplings of qnarks and 
charged leptons to the Higgs field are the only sources of flavour structure. The singular 
values of the Yukawa matrices, together with the vacuum expectation value (VEV) of 
the Higgs field determine the fermion masses, and the relative orientation between the 
up- and down-qnark Ynkawa matrices results in the Cabibbo-Kobayashi-Maskawa (GKM) 
matrix, responsible for charged flavour transitions in weak interactions. In the quantnm- 
field theoretical formulation of the SM, the Ynkawa matrices enter as coupling parameters 
in the Lagrange density. In the following, we will focus on the quark sector, where one has 

-/:Yukawa = YH QlH Uj^ + Q^H Dj^ + h.C. 
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( 1 . 1 ) 




Here and in the following Yjj and Yd denote the Yukawa matrices for up- and down- 
type quarks, Ql^Ur^Dr are the left-handed quark doublet and right-handed singlets, 
respectively, and H, H is the Higgs field and its SU{2) conjugate. The indices i, j = 1.. .Ug 
denote the quark generations/families (n^ = 3 in the SM). 

As all other couplings in the SM, after renormalization of ultraviolet divergencies, 
the Yukawa matrices in (1.1) are to be interpreted as effective parameters with the scale- 
dependence controlled by renormalization-group (RG) equations [1-4]. The structure of 
the RG equations and their solutions have been extensively studied in the past. In [5], the 
resulting one-loop RG evolution of the GKM matrix elements (in a given parametrization) 
has been studied, and approximate analytic solutions have been derived on the basis of 
the observed hierachies in quark masses and mixing angles in the SM. Generalizations 
to particular new physics (NP) frameworks have also been derived, notably for 2-Higgs- 
doublet models or supersymmetric extensions of the SM, see, for instance, [6, 7]. Recently, 
the effect of possible NP contributions has been studied in a model-independent way, by 
considering the RG effects from dimension-six operators in an effective field theory (SM- 
EFT) approach [8]. Finally, Bednyakov et al. [9] have recently computed the three-loop 
RG coefficients for the SM Yukawa matrices. 

The RG equations are usually formulated in matrix form, i.e. the scale-variation of the 
Yukawa matrices is given by a matrix polynomial of Yjj and Yr. Since the gauge sector 
of the SM is invariant under unitary field redefinitions for the individual quark multiplets, 
the RG equations have to transform covariantly under such changes of flavour basis (see 
below). This also implies a certain degree of redundancy in the RG equations, because from 
the 18 complex matrix entries in Yjj and Yd only 10 physical parameters are observables. 

In this paper, we will therefore reformulate the RG equations in terms of flavour 
invariants, i.e. objects constructed from Yjj and Yd which are independent of the choice 
of flavour basis. As has been shown in [10] from the basic algebraic principle of Hilbert 
series, one can define eleven polynomially independent flavour invariants for three quark 
generations. These fix the six quark masses, the three mixing angles and the sine and 
cosine of the CP-violating phase in the GKM matrix. As a corollary, using Cayley-Hamilton 
identities for matrix products (cf. [11]), this also implies that any flavour-covariant product 
of Yukawa matrices that appears on the right-hand side of the RG equations for Yjj or Yd 
can be reduced to a finite basis of flavour matrices with coefficients given as polynomials 
of flavour invariants. It is then a straightforward, though tedious, task to derive the RG 
equations for the set of flavour invariants. 

Although the RG equations for the flavour invariants contain the same information 
as the original flavour-covariant equations, the formulation in terms of flavour invariants, 
under certain circumstances, may be considered advantageous. For instance, the form of 
the RG equations is universal, not only for the SM, but also for all extensions that obey the 
principle of minimal flavour violation (MFV) in the technical sense of [12]. Furthermore, 
the hierachical pattern of masses and mixing directly translates into a well-defined power 
counting for (suitably chosen) flavour invariants, which can be exploited to simplify the 
RG equations. An attractive physical picture arises if one assumes these hierarchies to be 
associated to some dynamical NP mechanism that can be traced back to an effective poten- 
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tial which determines the flavour structures at low energies. Within the MFV framework, 
the potential itself will have to be formulated in terms of flavour invariants, and the mini¬ 
mization of the potential should generate VEVs for the flavour invariants that reflect the 
particular pattern of (sequential) flavour-symmetry breaking in the SM (see [13]). Recent 
studies along these lines can be found, for instance, in [14-17]. Finally, our approach could 
be extended and generalized to cases where there are additional flavour structures in some 
tensor representation of the SM flavour symmetry group. For example, these could show 
up as coupling constants in front of higher-dimensional operators in SM-EFT [18, 19], or 
as new spurion fields in the MFV framework [20, 21]. 

The remaining paper is organized as follows: In the next section, we will first discuss 
a toy scenario with only two generations (2G) of SM quarks. The simplifications in the 
2G case (no CP violation, closure of SU{2) matrices under multiplication, small number of 
polynomially independent invariants) allow us to introduce our approach in a very trans¬ 
parent way, perform all calculations analytically and illustrate the RG equations for the 
flavour invariants in a graphical way. To this end, we will first give convenient definitions 
for flavour invariants and basic flavour matrices. In terms of these, the general form of the 
RG equations for flavour invariants will be derived. We also present analytical and numer¬ 
ical solutions for the RG equations that can be obtained from exploiting SM-like flavour 
hierarchies in the one-loop approximation. In Section 3 we generalize our framework to 
the realistic case of three quark generations (3G). To keep the discussion transparent, we 
restrict ourselves to the one-loop approximation from the very beginning. Again, we derive 
the general form of the (one-loop) RG equations for the eleven flavour invariants, and dis¬ 
cuss their approximate solutions in the SM. We close the paper with a short summary and 
outlook in Section 4. Some technical details about the use of Gayley-Hamilton identities, 
the explicit form of the 3G flavour invariants, and the general form of the two-loop RG 
equations in the 3G case can be found in the appendices. 


2 Two Quark Generations 


As mentioned above, in this section, we restrict ourselves to two generations of left-handed 
quark doublets and right-handed up- and down-quark singlets in the SM. The gauge-kinetic 
terms of the SM Lagrangian are flavour-blind, and therefore independent unitary rotations 
of the quark multiplets define a flavour symmetry, 

^quark = U{2f/U{1)B ~ SU{2)q^ X U{2)u^ X U{2)d^ , (2.1) 

which is only broken by the Yukawa couplings in (1.1). Here we factored out a U{1)b sym¬ 
metry for baryon number conservation, which is unaffected by the SM Yukawa interactions 
[13]. (More precisely, we find it convenient to factor out a U{1)q^ transformation acting 
on the left-handed doublets only). In a particular flavour basis, the Yukawa matrices for 
up- and down-type quarks read 


Yu = 



Yd 


Ecabbibo 


f Vd 0 
\0 Vs 


( cos 6 sin 6 
— sin 9 cos 9 


f Vd 0 

V 0 2/5 



( 2 . 2 ) 
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(2.3) 


Under a change of flavour basis, the Yukawa matrices transform as 

>1/ ^ ^ Uq, Yd , 

where Vq^ G SU{2)q^ etc. 


2.1 Flavour Invariants 


In the 2G case, one can construct five flavour invariants that are polynomially independent 
(see e.g. [10] and references therein for the mathematical background). In the following, 
to set the stage for the 3G case to be discussed in Sec. 3, we will discuss the construction 
and properties of these flavour invariants step by step. 

From the Yukawa matrices Yjj and Yd, one can construct flavour invariants in terms 
of traces or determinants of matrix products constructed from the non-negative hermitian 
matrices 

U = YuY^ , D = YdY^ . (2.4) 

These transform as Vq^ U Vq^ and Vq^ D Vq^ under basis tranformations for the left- 
handed quark doublets. A convenient choice for non-negative invariants is 

h = tr[t/] =yl + yl>^, tr[T>] =yl + yl>f), 

h = det[[/] = ylyl >0 h = det[T)] = ylyl > 0, (2.5) 

and 

h = tr [U D] = {ylyl + ylyl) cos^ 6 + {ylyl + ylyl) sin^ 0 > 0 . (2.6) 


Apart from discrete ambiguties related to renaming the original quark fields in the flavour 
eigenbasis, they determine the four eigenvalues for the Yukawa couplings and the Cabibbo 
mixing angle. Invariants built from traces of higher powers of U and D are related to the 
above via Gayley-Hamilton identities (see appendix A and e.g. the discussion in [11]). Still, 
for the following discussion, we further define the polynomially dependent invariants 


13 = ^-213 = ^{yl - y, 


2^2 

u) 


_ -'2 


1 


I,^:^-2h=^-{yl-yl)\ 


(2.7) 


and 


h^h-^ = ^{yl- yl){yl - yl) cos{29 ), 


1 


h = hh - h - - Uu 


2 \ 2 / 


ids - Vd) sin {26). 


( 2 . 8 ) 


Triplet Matrices and Triplet Invariants: It is further convenient to divide the ma¬ 
trices U and D into singlet and triplet components with respect to the flavour group factor 
SU{ 2 )q^, 


U=^tv[U]l + U 3 , 


D = 1 triD] 1+ D 3 , 


(2.9) 


A third independent triplet matrix can be defined as 


A3 = -[U3,D3] 


( 2 . 10 ) 
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Polynomial Basis: For generic Yukawa entries, any 2x2 matrix M that transforms as 
Vq^ MVq^ under SU{2)q^ can be written as a finite polynomial of matrices from the set 
{1, C/ 3 , Ds, A 3 }. For those matrices, the following multiplication tables for symmetric and 
anti-symmetric products of matrices holds. 


{, } 

Us 

Ds 

As 

Us 

h 1 

h 1 

0 

Ds 

h 1 

/4I 

0 

^3 

0 

0 

1 


^[, ] 

Us 

Ds 

As 

Us 

0 

2A3 

-2Ds 

Ds 

-2A3 

0 

2Us 

As 

2Ds 

-2Us 

0 


This explicitly shows that the set {1, C/ 3 , 1 / 3 , A 3 } closes under matrix multiplication with 
prefactors that are polynomials of the flavour invariants . Here we defined the dual matrices 

C73 = |[Z/3,^3] = §t/3-|l?3, ^3 = |[A3,C/3] = |l/3-|t/3, ( 2 . 11 ) 

and 

As='-[Us,D3] = A3, ( 2 . 12 ) 

which can be obtained from the inverse of the metric 


as 


Gxy = tr[Y3Y3] 


(hh 0 ^ 

h h 0 

VO 0 1 / 


{X,Y = U,D,A), 


^ Us \ 


^Us^ 

Ds 

~ 2^ 

Ds 

\As) 


[as) 


(2.13) 


(2.14) 


From this we can read off the orthogonality relations between triplet matrices and their 
dual, 


tr [X 3 i 3 ] = ^ hxv (Y, Y = U,D,A). 


(2.15) 


This can be used, for instance, to decompose a generic SU{2)q^ triplet matrix as 


Ms 


2 tr[Y 3 M 3 ] 

2 ^ - 7 - ^3 

X=U,D,A -'5 


(2.16) 


Similarly, a matrix Mjj that transforms as a bi-doublet under SU{2)qj^ x SU{2)u^ can be 
decomposed as 


Mu 


tT[MuYjJ^] ^ ^ ^ 2_i4XsMuY^^^ 

^ X=U,D,A -^5 


(2.17) 
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and analogously for Un -H- D/j. Higher tensor representations of the flavour symmetry 
group and their expansion can be constructed from Notice that for generic matrices 

M 3 ^u,d, the coefficients in these expansions are enhanced by and respectively. 

In contrast, the MFV hypothesis assumes these coefficients to be of order 1 or smaller (see 
again [ 11 ]). 


2.2 Renormalization-Group Equations 
2.2.1 General Form 


The Yukawa matrices are subject to renormalization-group (RG) evolution. The generic 
form for the RG-running of the Yukawa matrices Yjj^d can be written in manifestly flavour- 
symmetric form (see e.g. [2]). Using the generic decomposition into basis matrices as 
discussed above, we thus write 


dYuif,) 

dlnfi 

clYpifi) 

dlnfi 


{ao{Ii,fi) 1 + ai(/i,/i) Us + D 3 + ia 3 {Ii,^i) A 3 )Yu{ij) , 

{bo{Ii,ld) 1 + bi{Ii,fi) D 3 + b 2 {Ii,iJ.) U 3 -ib 3 {Ii,fi) H 3 ) YD(fi). 


(2.18) 


Each of the coefficients Oj, bi depends on flavour invariants which arise from loop diagrams 
including additional Higgs-Yukawa couplings. (At one-loop accuracy, only terms at most 
quadratic in the Yukawa couplings can appear within the round brackets etc.) In the SM 
(or, in general, in constrained MFV models without additional sources of GP violation), the 
coefficients will be real polynomials of the flavour invariants.^ This immediately translates 
into RG equations for the matrices U and D, an from this we obtain 


dh 

dlnn 

dh 

dlnfi 


tr 


tr 


dU 

dlnfi 

dD 

dlnfi 


2 (oo h + aih + 02 h) , 
2 {bo h + bil4 + 62 h) ■ 


(2.19) 

( 2 . 20 ) 


In a similar way, one obtains the RG equations for the remaining invariants in a 
straightforward manner. The RG equations for the invariants I 3 , I 4 take a particularly 
simple form 


dls 

dlnfi 


4ao h , 


dh 

dlnfi 


460 / 4 . 


For the invariant h, we obtain 

= (2ao + 26o -I- aih -|- bih) h + 02 / 1/4 + 62 / 2/3 + (03 + 63 ) h 1 
and for the invariant h, we get 


( 2 . 21 ) 


( 2 . 22 ) 


— (4ao + 46o -I- 2aili -|- 261/2 — 2(03 -|- 63 ) h) h , (2.23) 

dlnfi 


^Furthermore, if weak isospin-violating corrections are neglected, the coefficients ai and bi will be related, 
see e.g. [3, 22]. 
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Figure 1. Illustration of the “phase space” (the region below the shaded area) for the normalized 


invariants x, y, z defined in the text. 

Discussion: From (2.21) and (2.23) we observe that some limiting cases in the phase 
space of flavour invariants are stable under RG evolution: 

• The case I 3 = 0: 

In terms of physical parameters, this corresponds to one vanishing eigenvalue in the 
up-quark sector, = 0 , and otherwise generic values for yd,s,c and 6 . 

• The case I 4 = 0: 

This corresponds to one vanishing eigenvalue in the down-quark sector, = 0, and 
otherwise generic values for yu,s,c and 9. 

• The case I 5 = 0; 

This corresponds to sin20 = 0, i.e. no mixing and otherwise generic Yjj and Yd; or 
degenerate eigenvalues in the up-quark sector (y^ = y^) or in the down-quark sector 
{Vd = ys), respectively. 

For illustration, we thus define normalized invariants (for Ii ^2 7 ^ 0), 


4 G Ah Ah 


(2.24) 



which take values in the unit interval [ 0 , 1 ], with the additional constraints 


(2.25) 


z < xy. 


This is illustrated in Fig. 1. As a consequence of the above observations, there will be now 
RG flow from the “phase-space” edges, defined by x = 0, y = 0, or z = 0, into the bulk. 
This can be understood as a consequence of a residual flavour symmetry. In contrast, the 
case z = xy is not protected by symmetry. A more detailed discussion of the residual 
flavour symmetries associated with this situation will be given in [23] (see also [24]). 
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2.2.2 Exploiting Flavour Hierarchies 


The RG equations simplify when one exploits flavour hierarchies in the Yukawa matrices. 
For instance, in a SM-like scenario, we can consider the limit where all but one Yukawa 
coupling, say yc in the 2G toy case, are small. In this case, the basis of triplet matrices in 
(2.18) can be reduced to t/3, and consequently only the coefficients ao, bo, ai, 62 are relevant 
to first approximation.^ The RG equations for the invariants in this approximation read 
(also using Is ~ /^/ 2 ) 


dh 
din y 

dh 

din y 

dh 

din y 


— (2ao + aih) h , 


— “^boh + 262-I5 ) 


dh 

din 

dh 


= 400/3 


din /i 

~ (4ao + 46o + 2ai/i) h ■ 


— 460/4 


(2.26) 


Solving for the four coefficients, leaves one general relation between the five invariants and 
their derivatives which can be written as 


This implies 


dh „ dh dh If ^ N 

^^2 —+ — (for yc > 2/M,s)- 
h -'i h 


hM /f(yo)/4(/xo) ^ 


(2.27) 


(2.28) 


or, in terms of Yukawa eigenvalues and the Gabibbo angle, ^ 


{ys/Vd - Vd/Vsf sin^(26') ~ const. (for y^ > y„,d,s)- (2.29) 


Putting in experimental values for the quark-mass ratio and the Gabibbo angle, the con¬ 
stant on the r.h.s. ranges between 60 and 90. 


2.2.3 One-loop Solutions in the SM 


To illustrate the numerical effect of the RG equations, we consider the one-loop RG coeffi¬ 
cients in the SM. The system of RG equations further simplifies if we neglect electroweak 
corrections, leading to the values summarized in Table 1. For the starting values of the 
invariants in the 2G case, we consider a toy model where we neglect the first generation 
in the SM, such that the large Yukawa couplings from the third generation lead to non¬ 
trivial effects on the r.h.s. of the RG equations. Exploiting again the hierarchies in the SM 
Yukawa entries, we then find 


^ ~ JL /2 _ n 

d In y ^ tt 


(2.30) 


complementary approach would perform the limit j/c yu,d,s from the very beginning and consider 
invariants under the reduced flavour symmetry only, see [25, 26]. 

®In models with texture zeros one typically relates the Gabibbo angle to the square root of yd/ys, see 
e.g. [27]. Therefore, such relations - in general - are not scale invariant in the limit of hierarchical Yukawa 
couplings. 
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«0 

bo 

ai = bi 

02 = &2 

8^2 

[h + h + \-:;) 

_ 

TT 

8^2 


_ 

TT 

3 

1677^ 

_ ^ 

1677^ 


Table 1. Non-vanishing RG coefficients in the SM as defined in the text, following from [5]. 
Here Ug denotes the number of generations, and electroweak contributions have been neglected for 
simplicity. 


Using the one-loop expression for the QCD /3-function, 

da. 


d In/i 


s 2 

— —— a 


/3o 

— a. , 
2 tt ^ ’ 


one obtains the explicit solution 

h ( m ) ^ h (no) G in, Ho), 

where we defined the RG-evolution function 


G{h, Ho) ■= exp 


9 dn , . 


'Mo 


V = 


9 




asin) 

asiHo) 


‘fj — 00 


-1 


dvr (/3o - 8 ) asin) 


(2.31) 


(2.32) 


(2.33) 


This coincides with [5], where the approximate RG flow of the top Yukawa coupling has 
been derived (with Ii{h) — (ytin))'^ pO) in [5] is defined as in our convention.) 

For the remaining invariants, using Uc ^ Vs ^ yu,d and | sin20| <C 1, we have 


dh 3 Aus 


dlnH Att^ \Ii 


(2.34) 


and 


d /3 15 

dluH 
d /4 9 

d In ^ Svr^ 
d/s _ 27 
dluH 

We see that once the RG-solution for R (p) has been constructed, the RG equations for the 
remaining invariants can be easily solved by separation of variables. Using the RG function 
G{h,Ho) defined in (2.33), we have 


/1/3- 

8 as - 

-/s 

TT 


dx 

3 , 


dln/r 


/ 1 / 4 - 

8 as y 

- 14 

TT 


dy 

rs-' 

dln/r 

3 , 


1J5- 

Idas 7 

TT 


dz 

rs-' 

dln/r 


(2.35) 


l2iH)^H^^^°h{Ho)[G{H,Ho)f\ 


(2.36) 


and 


x{h) - [G{h, Ho)] ^^^x{ho), yin) ^ [G{h, Ho)]^^^ y{Ho), z{h) ^ [G{h, Ho)]^^^ z{ho) ■ 

(2.37) 
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2.2.4 Numerical Illustration 


In Fig. 2 we provide illustrations for the one-loop RG flow of the combinations of flavour 
invariants x, y, z in the SM, and compare the exact numerical solutions with the approxi¬ 
mation in (2.37). We observe that — for the chosen numerical starting values — even for 
values as large as t = Infi/fiQ = 15, the differences between the exact and approximate 
solutions are always below 5%. 


1.04 

3 1.02 
N 

1.00 
"S' 0.98 
0.96 

0 5 10 15 0 5 10 15 



1.04 



1.02 



rsi 

-1.00 





"^0.98 



0.96 



t-to 


t-to 


Figure 2. Comparison of numerical [dots] and approximate analytical (2.37) [solid line] solution for 
combinations of flavour invariants, x{t)z{t) and y{t)/z{t), normalized to the values at t = Iny /> 
to = 0. The following starting values have been used: /i(to) = 1-0, ^ 2 (^ 0 ) = 0.1, x(to) = 0.005, 
yito) = 0.048, z(to) = 0.105. 


In Fig. 3 we illustrate the RG flow for the boundaries of the “phase-space” of flavour 
invariants, defined byz = 0,y = 0,x = 0, respectively, see the discussion in Section 2.2.1. 
Again, we have chosen a hierarchical scenario with hifJ-o) ^ h{fJ-o)- We observe that 


• The relation y{t)/z{t) ~ const, holds on the whole plane x = 0, which is in line with 
our derivation of (2.28) which only required Ii l 2 - 


• In contrast, x{t)y{t) ~ const, only holds in the vicinity of x rNj y r\j 0 (where yu 'C 
yc and yy <C ys) and for 9 near zero (which requires the solution with I^{yo) = 

— \J 13/4 — / 5 (/ro )-^5 shown on the left-hand side). 

• The same is true for x{t)z{t) ~ const.. 
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y(t) y(t) 




z(t) 



z(t) 





Figure 3. Numerical illustration of the RG flow at the “phase-space boundaries” for flavour 
invariants in the SM (2G, one-loop accuracy, neglecting electroweak gauge couplings). Top: z = 0; 
center: y = 0; bottom a: = 0. Each arrow indicates the RG flow from fg = 0 to t = 5; the starting 
values are again chosen as Ii(/ro) = 1.0, = O-l? ctsif^o) = 0.2. The plots on the left (right) 

are generated assuming / 5 (/io) = Tl/ I 3 I 4 — hifJ-o)- 
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3 Three Quark Generations 

For three quark generations in the SM, the flavour symmetry group to consider now is 

^quark = U{2>f/U{l)B ~ SU{ 2 ,)q^ X X U{ 2 ,)d^ . (3.1) 

The corresponding Yukawa matrices again transform as bi-doublets under a change of 
flavour basis, 

Yu ^ VQ.YuVi^ , Y,, ^ Vq^YdVI^ . 

In a particular flavour basis, they are given by 



0 

0 


0 

0 

Yu = 

0 Vc 0 

, Yo = Vckm 

0 

0 


0 

0 


\ 0 0 yb J 


In the subsequent analysis, it turns out to be more convenient to discuss the flavour invari¬ 
ants as a function of the CKM elements Vij without choosing a particular parametrization 
in terms of mixing angles which would directly reflect the unitarity of the CKM matrix 
VcKM- 

3.1 Flavour Invariants 

As discussed in [10], the SM quark sector in the 3G case can be described in terms of 10-|-1 
polynomially independent invariants, which determine 6 Yukawa eigenvalues, 3 mixing 
angles and the cosine and sine of the CP-violating phase (in a given parametrization of the 
CKM matrix). With a similar procedure as in the 2G case, we will now explicitly construct 
a convenient set for these 11 invariants from the non-negative hermitian matrices, 

U = YuY}j , D = YdYI , 

which now transform under the SU{3)q^ flavour symmetry. For later use we also define 
the adjoint matrices, satisfying 


U adj U = det U etc. 

With this, we can easily construct a complete set of polynomially independent positive 
semi-definite invariants. For the unmixed invariants, we define 

Ii = tr ([/) > 0 , I 2 = tr {D) > 0 , 

Is = tr (adj U) > 0 , I 4 = tr (adj D) > 0 , 

Iq = det (U) > 0 , Is = det (D) > 0 , . (3.3) 

These determine the six singular values of the Yukawa matrices. The GKM elements are 

then determined by mixed invariants which we define in a similar way. CP-even invariants 
can be chosen as 

Is = tr (C/H) > 0 , 17 = tr {D adj U) >0 , I 9 = tr ([/ adj H) > 0 , (3.4) 
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and 


/lo = tr (adj (C/L>)) > 0, (3.5) 

According to the discussion in [10] , there is an eleventh, CP-odd, invariant that cannot be 
expressed as a polynomial of the other ten invariants, as defined above. It is related to the 
Jarlskog determinant [28] and can be chosen as 

/fi = --det[[/,D] . (3.6) 

Explicit expression in terms of Yukawa couplings and CKM elements can be found in 
Appendix B. 


Octet Matrices and Octet Invariants: As in the 2G case, we can also construct 
basic flavour matrices as octet representations of the flavour group factor SU{3)q^. First, 
there are two polynomially independent octet matrices that are quadratic in the Yukawas, 
namely the traceless part of the matrices U and D (defined analogously to the 2G case), 

[/g = [/ _ 1 tr[17] 1, Ds = U-^tT[D]l. (3.7) 

For terms quartic in the Yukawas, we may define the octet part of the adjoint matrices, 

adj Us = adj U — - tr[adj U] 1 , adj Ds = adj D — - tr[adj D] 1, 
o o 

together with 

Ss^\{U,D}-^tT[UD]l, As='-[U,D]. (3.8) 

Similarly, we define 


Au^'- [adj U,D], [adj D, U] . (3.9) 

For generic Yukawa entries, the eight hermitian matrices as defined above provide a 
basis for octet matrices in S'C/(3 )q^. The symmetric (but non-orthogonal) metric defined by 
the traces of matrix products contains flavour invariants for the 3x3 case. It is summarized 
in Table 2. Here, the unhatted invariants are related to the hatted ones via 


h = tr[[/|] = hi - 2% , 


h = tr[C/8 adj Us] = 31^ 


hh 

3 


I4 — ti'[ 7 l|] — -I2 — 2/4 , 

Is = ir[Ds adj Ds] = 3/8 - , (3.10) 


and 


h — tr[C/8.D8] — I 5 -^ 


and 


/y = ti[Ds adj Us] = /y - 


hh 


h = tj:[Us adj Ds] = h- 


hh 


(3.11) 


(3.12) 
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tr[ ] 

Us 

Ds 

adj Us 

adjDs 

^8 

^8 

Au 

Ad 

Us 

h 

h 

h 

h 

h 

0 

0 

0 

Ds 

+ 

h 

h 

h 

h 

0 

0 

0 

adj Us 

+ 

+ 

Pi 

P 2 

P 3 

0 

0 

0 

adjOs 

+ 

+ 

+ 

Pa 

P 5 

0 

0 

0 

Ss 

+ 

+ 

+ 

+ 

ho 

0 

0 

0 

^8 

0 

0 

0 

0 

0 

ho 

Qi 

Q 2 

Au 

0 

0 

0 

0 

0 

+ 

Qs 

Qa 

Ad 

0 

0 

0 

0 

0 

+ 

+ 

Qb 


Table 2. Traces of basic octet matrices constructed from YjjY^ and Y^Yj^ as defined in the text. 


and 


If — tr[C/858] — I 7 — hh + 
I9 = tr[L)8'S'8] = h — hh + 


2/1/5 
3 ' 

2/2/5 
3 


Finally, one has 


T + rc c 1 hh — ho {h — hh) h hh + hh 
ho = tr[S’858j =- - -+- - -H-^-. 

Further polynomially dependent invariants that appear in Table 2 are given by 

7 / 


ho — tr[^8] — -^10 + 2/10 -^ 


(3.13) 


(3.14) 


(3.15) 


and 


21 . 




h 

4 


21 : 


Pi = tr[adj Us adj Us] = -2hh + ^ , Pa = tr[adj Ds adj Ds] = -2hh + ^ , 


and 


P 2 = tr[adj Us adj Ds] = ho - 


hh 


(3.16) 


(3.17) 
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and 


P 3 = tr[adj Us Ss] 



Furthermore, 


P 5 = tr[adj Ds S's] 



Qi = trf^s Au] 
Q 2 = trf^s Ad] 


h (-^10 ~ hl 4 ) 3/4/6 _ hh _ hh 
2 4 2 ^ ’ 

h (/lo — hlj) _ ‘^hh hh hh 

' 9 4 9 9 ’ 


and 


Qs = iT[Au Au] 
Qs = tr[^£) Ad] 


h {ho — hh) h ( 3/9 — / 1/4 + hh) {hh — h) h 
2 2 ' 2 ' 
Ihho-hh) Is{2>h-hh + hh) {hh-h)h 

9 “T cy “T ^ 


and 


— tr[^;7 Ad] 


{h-hh){ho-hh) 

4 ^ 

(/e - hh){h - hh) 

4 


(/? ~ hh){h — hh) 

4 

- 2 / 6/3 • 


(3.18) 


(3.19) 


(3.20) 


(3.21) 


In order to project onto the eight basis matrices one needs the inverse of the metric in 
Table 2. The explicit result is rather lengthy and not very instructive, and we therefore 
refrain from quoting it here. We checked however that the metric is not singular for generic 
Yukawa entries. 


3.2 One-Loop RG equations 

In the 3G case, again, any flavour matrix that arises as a flavour-covariant product of SM 
Yukawa matrices Yu and Yd can be written as a linear combination of a finite set of basic 
matrices (constructed from Yu and Yd) with coefficients given as polynomials of a finite 
number of flavour invariants (as a corollary to the discussion in [10]). The most general 
form of the RG equations then can be written as 


dYu 

dlnp, 

dYp 

dln/r 


(ao 1 + aiUs + a2 Ds + as adj Us + 04 adj Ds -I- 05 5*8 -I- iae As 
+iaY Au + iasAp) Yu{fi) , 

{bo 1 + biDs + b2Us + 63 adj Ds -/ b^ adj Us Yb^Ss- ibo As 
—iby Ad — ibsAu) Yd{p) , 


(3.22) 


As compared to the 2G case, the expressions for the RG equations of the flavour invariants 
derived from this general parametrization become rather lengthy. (The explicit structure 
of the two-loop expressions can be found in (C.l) in the appendix.) 
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3.2.1 General Form 

If we restrict ourselves to the RG equations at one-loop accuracy, we can write 

dYu 


dln/r 

dYp 

dlnfj, 


— (oq 1 -h l/g Y a2 Dg Yp^p ,), 

= (6o 1 + Dg + b2Ug + ...) Yp{fj,), 


(3.23) 


where the coefficients ao,^o are hrst-order polynomials of flavour invariants, and ai^ 2 )fei ,2 
are constant, see again Table 1. The RG equations for the quadratic invariants then take 
the same form as in the 2G-case, 

~ 2ao h + 2oi Ig + 202 I 5 j — 2&o I 2 + 26i I 4 -|- 262 h ■ (3.24) 

We remind the reader of the difference between the hatted and unhatted invariants, as de¬ 
fined in Section 3.1. For the remaining unmixed invariants, we also find simple expressions 


dig 

~ 4ao Ig — 2 ai Iq — 2 a 2 I 7 ■ 
at 


~ 46 o h — 2&1 Ig — 262 19 , 
at 


and 


dh (, f dig - 

— Doo Iq , -j- — 660 h ■ 

at at 


(3.25) 


(3.26) 


Notice that the last two relations — with our convention in (3.22) where the coefficients 
always multiply traceless matrices — are exact. The one-loop RG equations for 
the mixed invariants are determined as 


— ( 2 ao + 260) Ib + ( 2 ai -|- 262) h + (202 -I- 2bi) Ig . 


(3.27) 


together with 


dir -~~ 

= ( 4 ao -I- 26 o) Ij + ( 2 oi — 262) 


— (202 — 25i) /lo — Igl^ -|- 


— ( 46 o + 2 oo) Ig -I- ( 26 i — 202) 


I hh 

V 3 

2 / 2 / 7 ^ 
( hig 


-hh 




(3.28) 


{' 


— (262 — 2 ai) /lo — Igli + 


\ ^ 

2hlg \ 

3 ) 


-hh 


(3.29) 


and 


dlio 

dt 


— (4ao + 46o) 7io + ( 2 oi -|- 262 ) 


hho 


— hh 1 + (2 o2 -|- 2 bi) 


42^10 


-hh 


The RG equations for the Jarlskog invariant is simple, and to one-loop accuracy reads 

dl 


‘11 


dt 


— (600 -|- 660 + 2 oi h Y 2 b\ I 2 ) III • 


(3.30) 
ds 

(3.31) 
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3.2.2 Exploiting Flavour Hierarchies 

The RG equations again simplify when one exploits flavour hierarchies in the SM Yukawa 
matrices, which are also applicable to MFV extensions of the SM. For concreteness, we 
relate the scaling of the quark Yukawa couplings to the Wolfenstein parameter A in the 
CKM matrix, as it can be realized in Froggatt-Nielsen models [29] (see also [13, 20]), 
assuming 


Vi2 ~ A, V23 ~ A^ , V13 ~ A^ 


(3.32) 


and 


yt ~ , yh^ >? 1 ys ~ A® , yu,d ~ A® . 

Defining e = A^, the individual invariants scale as (see Appendix B) 


h = viycVu ~ e 


h ^ y? ~ , h- ytVc ~ , 

h^yl^e^, 74 ~ ylyl ~ e®, h = ytytyj ~ e 


11 

2 „. 2„,2 . ,16 


(3.33) 


(3.34) 


and 


/5-y?y6|Gifel2 = /ii2 + 0(e"), 

ho - yhlylyllVud? = hh + 0 (e ^^), 


h - yhlvllhubl 
h - yhbyl\ytd\ 


2 ,8 

~ e , 


2 ,11 

~ e 


(3.35) 


and 


hi ~ e 


16 


The leading terms in the (one-loop) RG equations are then identified as^ 


and 


and 


dh 


dl^ (, 2ai 


dh 


462 


— 1 260 + 7i ) h 


^ ~ ( 4ao + ^ /i ) Is , 


dh 


262 


~ ( 46o + — 7i ) 74 


d{h — hh) 
dt 


— ( 2 ao -I- 2 bo ^ 7 i-^ 7 i ) {h — I1I2) 


(3.36) 


(3.37) 


(3.38) 


dh f 

la 


dh ( 2ai — 262 , 1 

— =(4a„ + 2i,„ +-^- h]Ir 

— =. (4l,„ + 2o„ +-j 


(3.39) 


^ We do not include the power-counting for the gauge-coupling constants here as has been advocated in 
[ 22 ]. 
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and 

- (4«(, + 46„ + /.) (/.„ - m . (3.40) 

together with 

~ ( 6 ao + 660 + 2oi/i). (3.41) 

As in the 2G example, only the coefficients oq, bo, ai, 62 in (3.23) are needed in this approx¬ 
imation. Solving for the latter, one obtains 


oo dt 
60 dt 


6/1 3/3 ’ 

(i/2 (i/4 


oi dt ~ 


1 



62 dt ~ 


1 




(3.42) 


This leaves 7 relations that can be used to identify RG-invariant combinations of flavour 
invariants. 


and 


and 


and 


die , dll dis dis d/2 „ dh 

^ + — = 2 ^, —~ 2 ^, 
le Jg Jg 72 /4 


d{h - hh) ^ dh ^ dh ^ dh 

/5-/1/2 h ~ h /4 ’ 


dir , d/2 d/3 d/4 d/g d/2 d/i dh 

^ + — ^+2— -K2^, 

h ^2 /3 /4 /g 7 2 7l /4 


d(/io — -^ 3 / 4 ) d /3 dh 

— - — ~ , 

ho — hh h h 


dill , d /2 d/i d /3 d /4 

+ — ^ — + — + 2 — 

Al ^2 /l /3 /4 


(3.43) 


(3.44) 


(3.45) 


(3.46) 


Here each of the invariants is to be read as a function of (R, h^h^h)- As in the 2G case, 
the relations can be easily integrated, resulting in 


and 


~ const. ~ ^ ~ , 4^^ ~ const. ~ ^ ~ , 

{h? vl {h? 


(3.47) 


hh . Vs 2 

~ const. ~ n ~ e . 


hih — hh) yb|f 4 bp 


(3.48) 
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and 


hh ^ Vs 

~ const. ~ 

^2-^7 2/f) \^ub\ 


hlj 

hh 


~ const. ~ 


yl\Vu[^ 


(3.49) 


and 


ho — hh . |2 hhh 

—-~ const. ~ I Vus\ ~ e j -^ 

hh hhi 


~ const. ~ 


^yl 




e. 


(3.50) 


This explicitly shows, how the known simplifications for the RG solutions of quark masses 
and mixing angles that arise in the limit of large top-quark Yukawa coupling (see also [30] ) 
can be translated to the set of flavour invariants in a straightforward manner. 


3.2.3 One-loop Solutions in the SM 

As in the 2G-case, we can derive explicit solutions to the RG equations, using the one- 
loop expressions for the coefficients in Table 1 and the approximations from the Yukawa 
hierarchies discussed in the previous paragraph. With our definitions of flavour invariants, 
the approximate RG equations for the invariants h ,2 and h, 4 : looks identical to the 2G case 
in (2.30,2.34,2.35). As a consequence, we can again express the running of the 11 invariants 
in terms of the RG function G(/x, hq) defined in (2.33) from the evolution of the leading 
invariant 

/i(^) ~r/®/^o/i(^o)G(Ai,/io), V = ■ (3.51) 

as[yo) 

Defining normalized invariants as before (using a slightly different notation), we have 

X2iy) = ^ [G{n, X2{fJ-o) , 

h{y) 

X 3 {y) = - [Giy, xsiyo), 

Xi{y) = ~ xhyo) ■ (3.52) 

and the remaining scaling relations follow from (3.47-3.50). In this way, we recover the 
results for the approximate RG running of CKM mixing angles as discussed in [5]. 


Comparison with Harrison et al. In a paper by Harrison et al. [31] it has been 
highlighted that, within the SM, the one-loop RG equations exhibit two combinations of 


d ( iv[UD] \ 

dt \(det[C/D])^/^J 

I (tr [{UD)-^] {det[UD]f/^] 


d h 


dt 

ho 

'dt ihls^h 


flow. 


= 0 (SM@l-loop), 

(3.53) 

= 0 (SM@l-loop). 

(3.54) 
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In our notation, we have 


d Is _ 2 (ai + 62) -I7 + 2 (02 + 61) Ig 

dt (leTsy/^ ~ {det[UD]y/^ 

d ho 2 (ai + 62) (^-^1-^10 — 3/4/6^ + 2(02 + 61) (^/2-fio — 3/3/8^ 
dt{%%Yh~ 3{det[UD]f/^ 


(3.55) 

(3.56) 


This indeed vanishes for ai = —62 = — 0.2 = bi which holds within the SM, see Table 1. 


4 Summary and Outlook 

From the experimental as well as form the theoretical side (see e.g. the reviews in [32-35]), 
the quark flavour physics program is currently entering the precision era. The goal is to hnd 
hints to physics beyond the Standard Model (SM) from dedicated experiments, notably 
LHCb and BELLE 11. Still, the answer to the flavour puzzle itself may reside at extremely 
high scales, possibly as high as the Planck scale. In any case, the determination of flavour 
observables occurs at low energies, and thus for any comparison with “new physics” models 
one needs to include the renormalization-group (RG) running of the flavour parameters 
in a given theoretical framework. In principle, there are various roads to discuss this. 
On the one hand, one can consider the entries of the 3x3 Yukawa matrices and study 
their RG evolution; but these depend on an arbitrary choice of basis in flavour space. 
Alternatively, one can use the physical parameters, i.e. the six quark masses together with 
four independent GKM parameters to describe quark mixing and CP violation in weak 
interactions; but these have rather complicated relations to the Yukawa couplings. 

In this paper, we have chosen an intermediate point of view and considered simple 
combinations of Yukawa couplings that are independent of the orientation of the flavour 
basis. In terms of these flavour invariants we have formulated RG equations which are basis 
independent and allow for a transparent implementation of flavour hierarchies as observed 
in the SM or its minimal-flavour-violating (MEV) extensions. Expanding systematically in 
small parameters, we have also constructed simple analytic solutions for the RG evolution 
of a set of polynomially independent flavour invariants. 

Discussing the RG flow in terms of flavour invariants may be advantageous to discuss 
models with dynamical flavour symmetry breaking, where the Yukawa couplings emerge 
as vacuum expectation values (VEVs) of some scalar flavon helds. The scalar potential 
generating these VEVs will be constructed in terms of polynomials of flavour invariants 
of a given canonical mass dimension. In MEV-like constructions (see e.g. [36]), these can 
be reduced to the set of invariants discussed in this work. More complicated situations 
arise if one implements the spontaneous breaking of a gauged flavour symmetry on the 
level of renormalizable interactions. This leads to an “inverted-MFV” scenario, where 
the fundamental flavour invariants are approximately given as polynomials of the inverse 
Yukawa matrices [37]. Even more complicated relations can arise in a recently proposed 
model with dynamical flavour-symmetry breaking with a unification scheme according to 
Pati and Salam [38]. While the general form of the RG equations (3.22) will remain the 
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same, the coefficients will have a more complicated dependence than in MFV scenarios. In 
any of these cases, the renormalization-group flow of the invariants is needed to constrain 
the theoretical NP parameters at a high scale from flavour observables at low scales, and 
eventually give us some clue on the solution of the flavour problem. 
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A Cayley-Hamilton Identities 
A.l Two-Generation Case 

The Cayley-Hamilton identity for 2x2 matrices M reads 

0 = M2-tr[M]M-hdetMl. (A.l) 

Taking the trace and solving for det M, one obtains 

detM = ^ (tr2[M] - tr[M2]) . (A.2) 

Multiplying (A.l) with M~^, and solving for adj M = M“^detM, one obtains 

adj M = tr[M] 1 — M => tr[adj M] = tr[M] . (A.3) 

Inserted back into (A.l) yields 

= tr[M] M - det[M] 1. (A.4) 

For traceless matrices, this further simplifies to 

^ tr[M2] 1 (tr[M] = 0 ). (A.5) 

Therefore any power of 2 x 2 matrices M can be reduced to the basis {1,M} with co¬ 
efficients built from polynomials of trM, trM^ which are invariant under unitary basis 
transformations. 

For matrices Y which transform under bi-unitary transformations, Eq. (A.3) general¬ 
izes to 

adj y = yt (tr[M] 1 - M) {M = yyt). 


(A.6) 




A.2 Three-Generation Case 


The Cayley-Hamilton identity for 3x3 matrices M reads 

0 = - tr[M] ^ (tr2[M] - tv[M^]) M - detMl. 

Taking the trace and solving for det M, one obtains 


1 


detM = - I tr[M'^]-tr[M]tr[M^] H—tr'^[M] 


(A.7) 


(A.8) 


Multiplying (A.7) with M and solving for adj M = M ^ det M, one obtains 

adj M = — tr[M] M + ^ (tr^[M] — tr[M^]) 1 

=> tr[adj M] = ^ (tr^[M] — tr[M^]) . (A.9) 

Inserted back into (A.7) yields 

= tr[M] - tr[adj M] M + det[M] 1. (A.IO) 

For traceless matrices, this further simplifies to 

= - trfM^] M + - tr[M3] 1 (tr[M] = 0). (A.ll) 

2 3 

Therefore any power of 3 x 3 matrices M can be reduced to the basis {1,M, M^} with 
coefficients built from invariants that are polynomials of trM, trM^, trM^. 

Similarly as before, for matrices Y which transform under bi-unitary transformations, 
Eq. (A.9) generalizes to 

adj y = (tr[adj M] 1 - tr[M] M + M^) (M = yyt). (A.12) 


B 3G Flavour Invariants, Yukawa Couplings and CKM Elements 

For our convention to define 10+1 polynomially independent flavour invariants, the explicit 
expressions in terms of Yukawa couplings and mixing angles read as follows. The quadratic 
invariants are 


/i=tr[t/]= ^ l2^t4D]= yl (B-1) 

i=u,c,t j=d,s,b 

Again, Ii and I 2 quantify the overall size of flavour-symmetry breaking in the up- and 
down-quark sector, respectively. Quartic invariants appear as 

I 3 = tr[adj U]= Vi ^ = tr[adj D] = ^ ^ , (B.2) 

i=u,c,t j=d,s,b 
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and 


h = ir[uD] = Y. H y^yj ’ (B-3) 

i=u,c,t j=d,s,b 

where we have defined = VcUt, Vd ~ ylvt Continuing with the sixth-order invariants, 
we have 

% = tr[C/ adj [/] = 3 det [/ = 3 ylylvt , h = tr[L» adj D] = 3 det = 3 y^y^yl , (B.4) 
and 


/7 = tr[adj [/£>] = Y Y 

i=u,c,tj=d,s,b 

% = tic[U adj D]= Y Y y'i^i 


The eight-order invariant reads 


/lo = tr[adj [/adj D] = Y Y y‘iy‘j\^ii\^ 

i=u,c,t j=d,s,b 


(B.5) 


(B. 6 ) 


Finally, the CP-odd invariant 


hi = = - {yt - yl){yt - yl){yl - yl){yl - yl){yl - yl){yt - yl) Im > 

(B.7) 


is proportional to the Jarlskog determinant [28]. 


C Two-Loop RG equations for 3G Flavour Invariants 


The two-loop approximation for RG equations of the quark Yukawa matrices Yjj and Yd 
in (3.22) is obtained by keeping factors that are at most quartic in the Yukawa couplings, 
i.e. neglecting the contributions with the flavour matrices Ajj^d, 


dYu 

dln/i 

dYp 

dln/r 


(ao 1 ai C/s 02 T >8 + asadj C/g -h 04 adj Dg + a^S% + ia% ^g ...) Yu{yL) , 

(Co 1 + Cl Dg -|- 62 C/g -|- C 3 adj I/g -|- C 4 adj C7g -|- C 5 5g — ib% Tg -|-...) Y^iy) ■ 

(C.l) 


From this ansatz, it is straightforward - though tedious - to calculate the two-loop RG 
equations for the eleven flavour invariants. First, we have 


—r- ~ 2ao Ii + 2ai -|- 2a2 I5 + 2a^ Iq + 204 Ig Y 2 a^ Ij , 

at 


~ 2Co I2 + 2 bi Ii + 2C2 Is -|- 2C3 /g -|- 2C4 Ij + 2C5 Ig , 

at 


(C.2) 
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where here and in the following the abbreviations for the combinations of flavour invariants 
are the same as in Table 2. Then 


~ 4ao Is — 2 ai Jg — 2 a 2 I 7 — 2 as Pi — 204 P 2 — 2 as Ps , 

— 46 o Ia ~ 2bi Is — 262 Ig — 263 P4 — 264 P2 — 265 Ps (C- 3 ) 

and 

For the mixed invariants, one obtains 

~ (2ao + 260) I5 + (2ai + 262) Ij + (2a2 + 2bi) Ig 

+ (2a3 + 264) Ps + (2a4 + 2bs) P5 + (2a5 + 265) pg + (2a6 + 260) Po ; (C. 5 ) 


dt 

dh 

dt 


and 


dp 

— ( 4 ao + 2 bg) p + (2ai — 262) 


PP 


+ (203 — 264) ^2/5/6 — 3/5/6 ~ 


3 

‘II 3 I 7 

3 

hh 


— hh — (2a2 — 26i) pg — pp + 


2/2/7 \ 

3 J 


+ (2a4 — 263) 


hP 


-hh 


+ (205 — 265) (2/4 — h) h ~\ ——I" Qi + ( 2 ( 3.6 + 260) Qi, 


— (2ao + 46o) Ig — (2a2 — 26i) 


3 
hlg 


— (2a4 — 263) 2/5/3 — 3/5/3 — 


3 

2/4/9 


— Ph — (262 — 2ai) pg — hh + 


6 


2 / 1 / 9 ^ 


— (2a3 — 264) 


'hh 


-hh 


— {2as — 265) [ (2/3 — if) h + — Q2 ] — (206 + 260) Q2 , 


(C.6) 


and 

dpg 

dt 


— (4:00 + 46o) -^10 + ( 2 ai + 262 ) 


/i/io 


~ hh 1 + (2o2 + 261) 


hPo 


-hh 


(203 + 264) (/e {T9 - hh) + 


— (2o4 + 2 bs) {p — hh) h + 


2/4/10 \ 


+ (205 + 265) 


hPo 


— 3/0/3 — Q4 — (200 + 260) Q4 • 


(C.7) 
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and 


dl 


11 


dt 


~ ^6ao + 660 + 2ai/i + 262/2 ~ ~ 

+ i (05 + ^5) ^-(06 + be ) /s^ 


(C.8) 
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